Universality of ac-conduction in anisotropic disordered systems: An effective medium 

approximation study 
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Anisotropic disordered systems are studied in this work within the random barrier model. In 
such systems the transition probabilities in different directions have different probability density 
functions. The frequency-dependent conductivity at low temperatures is obtained using an effective 
medium approximation. It is shown that the isotropic universal ac-conduction law, u — a In cr, is 
recovered if properly scaled conductivity (a) and frequency (u) variables are used. 
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I. INTRODUCTION 

In the past years, ac-conduction in isotropic disordered 
systems has been extensively studied ,L2i£»£i£i2iiLiiiLii 
Several experiments support the existence of an univer- 
sal function governing the conductivity-frequency rela- 
tion in a variety of materials, in either electronio 5 *^ 7 " or 
ioniaaiSiSiiiiiSttii systems. This universality is also sup- 
ported by several theoretical studies^ One of the simplest 
models for studying conduction in disordered systems is 
the random barrier model, in which the energy barriers 
joining sites of a given network are selected at random 
from a given probability density function (PDF) . For this 
model the universality of ac-conduction at low temper- 
atures is well established. Dyre and coworkers studied 
the isotropic random barrier model within the effective 
medium approximation (EMA)m£&H They found that 
an universal function is arrived at in a low-temperature- 
small-frequency expansion; which corresponds to the so- 
lution of the equation 



= a In < 



(1) 



where a = er(u)/er(0) and u = u(3\n(3 / 'cr(O) are scaled 
conductivity and frequency variables respectively (with 
(3 oc (3 — 1/ksT, see below). This equation is univer- 
sal in the sense that it does not depend on the char- 
acteristic disorder of the system. The same equation 
had been previously arrived at with other approxima- 
tions, e.g. the macroscopic^ 2 ^ and hopping models^* In 
addition, other approaches such as the percolation path 
approximation, the diffusion cluster approximation and 
the velocity auto-correlation method, predict universal 
functions which present a better data collapse^' 7 How- 
ever, EMA offers a simple systematic tool to study ac- 
conductivity of disordered systems analytically which, in 
addition, gives the expected universal behavior. 

The conductivity properties of anisotropic disor- 
dered systems have also attracted attention in the last 
years, 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 Two relevant examples 
of anisotropic disordered systems are the superconduc- 
tor cuprates and porous reservoir rocks. In the first case, 
conductivity properties are strongly anisotropic, with a 
remarkable difference between the conductivity in the 



ab plane and along the c axis. In the second example, 
a relation between permeability and electrical conduc- 
tivity in isotropic fluid-saturated porous media is well 
establishedjS&iSSiS 7 . and a universal behavior for the dy- 
namical permeability, analog to Eq. QJ, was observed 
numerical and experimentally^^ 9 " Since anisotropy is a 
key characteristic of porous media and fractured rock^ 
the behavior of anisotropic frequency-dependent conduc- 
tivity in disordered media and its relation with the per- 
meability tensos^i is of key interest. 

In view of these and others examples of anisotropic 
disordered media, is necessary to find a relation be- 
tween anisotropic conductivity and frequency. The ap- 
proach used here is to use the ideas and concepts used 
in isotropic problems and apply them to anisotropic sys- 
tems. In this sense, the main purpose of this work is 
to extend and generalize the universal law, Eq. QJ, 
to anisotropic systems. An anisotropic random bar- 
rier model in two-dimensions is studied within an ef- 
fective medium theory. The analytical procedure used 
by Dyre* is closely followed in order to obtain a consis- 
tent low-temperature-small-frequency expansion for the 
anisotropic ac-conductivity. Here, from the analysis of 
the low temperature limit, the isotropic universal func- 
tion is recovered for properly scaled variables relating 
the conductivity in both directions. The paper is or- 
ganized as follows. In Sec. ^ the anisotropic random 
barrier model is described and previous zero-frequency 
results are summarized. Section IIIII briefly presents the 
main features of the anisotropic EMA and in Sec. IIVI 
the low-temperature-small-frequency expansion is per- 
formed. Finally, a summary is presented in Sec. 



II. BACKGROUND 

In the anisotropic random barrier model considered 
here, equal energy minima form a square lattice with 
its four nearest-neighbors separated by energy barriers, 
whose heights are randomly selected from anisotropic 
PDFs. Let 1 and 2 be the main directions of the 
square lattice, and the PDFs in each direction pi(Ei) 
and pi{E-i). Once the energy barrier E a p between two 



2 



nearest-neighbor sites a and (3 is selected from the cor- 
responding PDF, the transition rate between sites a and 
j3 is given by oj a p = u>oexp(— /3E a p), where uiq is the 
constant jump rate and (3 — 1/ksT is the inverse tem- 
perature, with ks being the Boltzmann constant. Since 
E a p = E/3 a , the forward (a — ► (3) and backward {[3 — > a) 
jumps have the same transition rate. The energy PDFs, 
pi(Ei) and p 2 (E 2 ), are related to the transition rate 
PDFs, v\(uj\) and ^2(^2)7 through the corresponding 
transformation of random variables. 

The zero- frequency conductivity of the anisotropic ran- 
dom barrier model in two dimensions was recently stud- 
ied by using an anisotropic generalization of the EMAiSi 
The low temperature conductivity in each direction was 
shown to follow Arrhcnius laws with the same activation 
energy E Cl which is determined by the anisotropic per- 
colation properties of the lattice. For the square lattice 
studied, the bond percolation threshold is the critical sur- 
face pi +P2 — 1 = 0, 31 " where pi represents the probability 
of having a conducting link between two nearest-neighbor 
sites in the i direction. This implies an activation energy 
given by24 

/ Pl (E 1 )dE 1 + p 2 {E 2 )dE 2 = l. (2) 
Jo Jo 

The zero- frequency conductivities, o~i(u — 0), in each 
direction are thus given by 



(72(0) = 7 2i^oa 2 e-^% 
where a is the lattice constant and the prefactor 



(3) 



E,. 



-1 
712 = 721 = 



/ P i{E 1 )dE 1 



E a 

J p 2 (E 2 )dE 2 




(4) 



Note that, at low temperatures, the anisotropic character 
of the system reflects only in the prefactors 712 and 721 
of the zero-frequency conductivity. 



III. ANISOTROPIC EMA 

The EMA consists in averaging the effects of disorder 
by defining an effective medium with effective transition 
rates, which depend on the Laplace variable u. These 
effective transition rates are self-consistently determined 
by the requirement that the difference between the prop- 
agator of the impurity and homogeneous problems should 
average to zeroi 32 i 33 i 34 i 3 ^ 3 ^ 3 Ti 3 ?i 3 ?i 4 9 In anisotropic prob- 
lems, two effective transition rates, one for each direction, 
are introduced. The effective frequency-dependent con- 
ductivities of the disordered system, <j\ (u) and a 2 (it) , are 
proportional to the effective transition ratesii^ii A ra- 
tionalized unit system is used where all the prefactors 



are absorbed in the definition of effective conductivities, 
so that they are equivalent to effective transition rates^ 
These effective conductivities are then determined by two 
self-consistent conditions ii^il 



(7i - LOi 



l + 2(a 1 -u 1 )[G 1 (u)-G°(u)} 



0. 



0"2 — Cl>2 \ n 

1 + 2 (a 2 - u*) [G*(u) - G°(u)] / Mua) = ' 



(5) 



Here, G 1 ^ and G° represent the non-perturbed 
anisotropic Greens functions related to the probabilities 
of moving from the origin to one of its nearest neighbors 
in the 1(2) direction and the return probability, respec- 
tively. The impure bond connects two nearest neighbor 
sites of the lattice whose transition rates are equal to u)\ 
if the impure bond lies in the 1 direction and u> 2 if the im- 
pure bond is in the other direction. The angular brackets 
denote averaging over the corresponding transition rate 
PDFs. 

The real and imaginary parts of the frequency- 
dependent conductivity are obtained considering that 
the Laplace frequency u is actually an imaginary fre- 



quency related to the real frequency by u 



4.38 



How- 



ever, the Laplace frequency picture has proved to be 
a useful simplification for studying frequency-dependent 
conductivity!^ 



IV. UNIVERSAL FREQUENCY-DEPENDENT 
CONDUCTIVITY 



In this section a low-temperature-small-frequency 
expansion 4 of the set of Eqs. (J5J is performed. In the 
following, the subscripts i,j = 1,2 are used to represent 
the two different directions of the lattice, noting that one 
of the self-consistent conditions in Eqs. (jHJ) is obtained 
interchanging the subscripts i and j in the other equa- 
tion. To the lowest order in u, the difference between 
the Greens functions in the i direction appearing in the 
denominator of Eqs. © may be written asi£ 



G u - G l 



fij , U 9ij 
~ 2(7, 2(7, ' 



with 



and 



9ij 



A-K^UiO-j U ((Ti + (Jj ) 



(G) 



(7) 



(8) 



Using this expansion for the Greens functions near u = 
0, the frequency-dependent self-consistent condition Eqs. 
(0 read 



{fij - UQij) 1 - 1 



0, (9) 
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In order to emphasize the energy dependence of the tran- 
sition rates, one may average over the energy distribu- 
tions and write, rearranging terms in the previous equa- 
tion, 



1 



fij ^Qij 



(fij - U9ij) 1 - 1 



(. 10 ) 

Two regimes may be identified for each direction: 

<7i, and uji (Ei) > 



mi) 

24 



1 



uj % {Ei) < I (f i:i 

(fij - ug l j)~ 1 - 1 
two cases, Ef(u), is defined by 

u>i[Ef(u)]= (fij-ugijY 1 -! 

which yields 

E a i (u) = ~\n{[(f i j-ug i j)- 1 -l 



aim* The energy separating these 



(11) 



UJ Q - 



-}• (12) 



Using these different ranges of energy, separated by E? , 
and the fact that transition rates and energies are re- 
lated through an Arrhenius law, the integral of the en- 
ergy average in Eq. 11U|) may be approximated, and the 
self-consistent condition, Eq. (JTJ3J, then reads 



Pi(Ei 



E?(u) 

or, equivalcntly. 



(fij ~ ugij) 1 - 1 



_ dE , = f.j ».'/,,, (13) 



J Pi(Ei)dEi = 1 - /, 



(14) 



EUu) 



In the zero-frequency case, Ef becomes the same for 
the two directions^ and is given by Eq. J2J) with Ef(u = 
0) = E c . Noting that for u = 0, fij is given by 



fij = J pi(Ei)dEi 



(15) 



and inserting Eq. (|15f) into Eq. (|14f> one obtains 

E c 

J pi(Ei)dEi = ugij, 

Ef(u) 

which may be approximated by 

e c 

Pi(Ei)dEi~ qi [E?(0)-E?(u)], 

E?(u) 



(16) 



(17) 



with q, = pi [E?(0)] = pi [E c ]. Evaluating Ef(u) from 
Eq. (O and combining Eq. (T5Jl with Eq. (JT^ one 
obtains, to the lowest order in u, 



ugij 



Si 
P 



In 



{fij ~ ug^y 



1 



(fr^- 1)^(0) 



m ■ 



p ai(oy 



(18) 



where <7i(0) is the zero-frequency conductivity in the i 
direction, and is given by Eq. J3J). Then, using the defini- 
tion of gij, Eq. |(SJ), the previous equation my be written 
as 



hi 



<7i(0) 



= Pi 



M(Ti(T, 

■ In ■ 



y/aiaj u(ai+aj)' 



(19) 



with /3j = f3/(Anqi). This last equation gives the gen- 
eral frequency-dependent conductivity for all tempera- 
tures and in the small frequency limit. 

In order to obtain the low temperature limit for the 
frequency-dependent conductivity, scaled conductivity 
and frequency variables for each direction are introduced, 
namely 



(20) 



\M(oK(o)' 



and 



Pi In Pi 



v*m*j¥) 

With this scaled variables Eq. (|19fl may be written as 



(21) 



ln< 



2 (n(0) 



\foiOj ln/3j 



In 



64cri(7j 
Ui(ai + Sj) 



+ In fa + In ( In /?, 



(22) 



By taking the low temperature limit p — > oo for fixed u~i 
and (Ti, the following set of coupled equations is arrived 
at for the scaled conductivities as functions of the scaled 
frequencies: 



ui = \J (Ti(T2(lncri -(- In 721 ) 
U2 = y/a^S\(hia^ + In 712) 



(23) 



where Eqs. © and (0} were used. 

This set of equations represent the complex relation 
between the conductivity in each direction. Although 
they may be regarded as universal equations for the 
scaled conductivities, the use of two different frequencies, 
one for each direction, is not suitable for a frequency- 
dependent description. In addition, the terms containing 
712 still depend on the specific PDFs used. Still, Eqs. 
(|23|l are useful to obtain the frequency-dependent con- 
ductivities: given the PDFs for each direction, E c and 
712 are calculated through Eqs. (0) and (@J, then through 
the definitions of Ui and Pi the scaled conductivities in 
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each direction may be calculated numerically by solving 
the set of coupled Eqs. lj2lfl) . 

As the analytical derivation of Eqs. (j23|l closely fol- 
lows the previous derivation of Eq. from EMAj4 the 
isotropic result is obviously recovered by setting the same 
PDF, p(E), for the two directions, which implies a\ — 
02 = o = cr(u)/cr(0) and u{ = U2 = u = u/3 In j3/cr(0), 
with j3 = (3/[4np(E c )]. However, by adding the two Eqs. 
(|23[) and setting u — (u± + u^) /2 and a — expres- 
sion Q is obtained again. This is a non-trivial result 
and establish that although the anisotropic conductivi- 
ties are coupled through the complex relation given by 
Eqs. I|23(l . the frequency dependence of the geometric 
mean conductivity, a, in the anisotropic problem is sim- 
ple given by Eq. Q). This result given by the EMA 
should be tested by more complex theories. 

V. SUMMARY 

In summary, by means of the frequency-dependent 
EMA, conductivity in an anisotropic random barrier 
model has been studied. It was shown that in a small- 
frequency expansion, the low temperature limit is char- 
acterized by an universal law, i.e. independent of the 
anisotropic PDFs, relating scaled conductivity and fre- 
quency variables. This scaled quantities are obtained by 
a proper combination of the conductivities and energy 
properties in each direction. Although the universal law 



is obtained for a two-dimensional system, it is expected 
to hold in three-dimensions for appropriately scaled vari- 
ables. Direct comparison with experiments on the con- 
ductivity of superconductor cuprates is not possible at 
present because the available data corresponds, to the 
best of the author's knowledge, to conductivity measure- 
ments in only one of the anisotropic directions of the 
system. Unfortunately, for a comparison with the scal- 
ing function, the conductivity in both directions should 
be available. The results presented can be of relevance 
for a complete theory of the anisotropic permeability of 
porous reservoir rocks. Given the relation between the 
anisotropic conductivity and permeability tensorsj2! the 
dynamical permeability in each relevant direction can be 
obtained through Eqs. (|23|l . Finally, is worth noting that 
a perfect agreement between the EMA universal law and 
experimental or simulation data should not be expected, 
as this is the case even for isotropic problems. 5-6 Other 
theoretical and simulation methods were shown to pre- 
dict a better universal law, which collapses experimental 
data from various disordered systems! 5 * 6 " However, EMA 
still provides a simple analytical tool for a first explo- 
ration of the properties of ac-conduction. In the present 
work, a first insight on the emergence of an universal law 
for anisotropic disorder systems has been presented. 
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